Theory of cosmological perturbations in an anisotropic universe 
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Abstract 

This article describes the theory of cosmological perturbations around a homogeneous and 
anisotropic universe of the Bianchi / type. Starting from a general parameterisation of the per- 
turbed spacetime a la Bardeen, a complete set of gauge invariant variables is constructed. Three 
physical degrees of freedom are identified and it is shown that, in the case where matter is de- 
scribed by a scalar field, they generalize the Mukhanov- Sasaki variables. In order to show that 
they are canonical variables, the action for the cosmological perturbations at second order is de- 
rived. Two major physical imprints of the primordial anisotropy are identified: (1) a scalar-tensor 
"see-saw" mechanism arising from the fact that scalar, vector and tensor modes do not decouple 
and (2) an explicit dependence of the statistical properties of the density perturbations and gravity 
waves on the wave-vector instead of its norm. This analysis extends, but also sheds some light 
on, the quantization procedure that was developed under the assumption of a Friedmann-Lemaitre 
background spacetime, and allows to investigate the robustness of the predictions of the standard 
inflationary scenario with respect to the hypothesis on the symmetries of the background space- 
time. These effects of a primordial anisotropy may be related to some anomalies of the cosmic 
microwave background anisotropies on large angular scales. 
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Inflation [j], 0] (see Ref. [3| for a recent review of its status and links with high en- 
ergy physics) is now a cornerstone of the standard cosmological model. Besides solving the 
standard problems of the big-bang model (homogeneity, horizon, isotropy, flatness,...), it 
provides a scenario for the origin of the large scale structure of the universe. In its sim- 
plest form, inflation has very definite predictions: the existence of adiabatic initial scalar 
perturbations and gravitational waves, both with Gaussian statistics and an almost scale 
invariant power spectrum [U, 0]. Other variants, which in general involve more fields, allow 
e.g. for isocurvature perturbations 0], non-Gaussianity jl], and modulated fluctuations 
All these features let us hope that future data will allow a better understanding of the details 
(and physics) of this primordial phase. 

The predictions of inflation are in agreement with most cosmological data and in par- 
ticular those of the cosmic microwave background (CMB) by the WMAP satellite |6|. The 
origin of the density perturbations is related to the amplification of vacuum quantum fluc- 
tuations of a scalar field during inflation. In particular, the identification of the degrees 



of freedom that should be quantized (known as the Mukhanov- Sasaki variables 10|), has 
been performed assuming a Friedmann-Lemaitre background spacetime [||. This means 
that homogeneity and isotropy (and even flatness) are in fact assumed from the start of 
the computation. In the standard lore, one assumes that inflation lasts long enough so that 
all classical inhomogeneities (mainly spatial curvature and shear) have decayed so that it 
is perfectly justified to start with a flat Friedmann-Lemaitre background spacetime when 
dealing with the computation of the primordial power spectra for the cosmologically ob- 
servable modes. This is backed up by the ideas of chaotic inflation and eternal inflation 0] . 
Note however that a (even small) deviation from flatness ll[ or isotropy [l2| may have an 
impact on the dynamics of inflation. It would however be more satisfactory to start from an 
arbitrary spacetime and understand (1) under which conditions it can be driven toward a 
Friedmann-Lemaitre spacetime during inflation and (2) what are the effects on the evolution 
and quantization of the perturbations. 

The first issue has been adressed by considering the onset of inflation in inhomogenous 
and spherically symmetric universes, both numerically in Ref. [l3| and semi-analytically in 
Ref. [14j. The isotropization of the universe was also investigated by considering the evolu- 
tion of four- dimensional Bianchi spacetimes (TBI . [l6l . 17j and even Bianchi braneworld [l8[ . 



No study has focused on the second issue, i.e. the evolution and the quantization of per- 
turbations during a non-Friedmannian inflationary stage, even though the quantization of 
test fields and particle production in anisotropic spacetime has been considered |19j . Such 
an analysis would shed some light on the specificity of the standard quantization procedure 
which assumes a flat Friedmannian background (see however Ref. |20j|). 

From an observational perspective, a debate concerning possible anomalies on large an- 
gular scales in the WMAP has recently driven a lot of activity. Among these anomalies, we 
count the lack of power in the lowest multipoles, the alignment of the lowest multipoles, and 
an asymmetry between the two hemispheres (see e.g. Refs. 2l|). The last two, which point 
toward a departure from the expected statistical isotropy of the CMB temperature field, 
appear much stronger. Various explanations for these anomalies, besides an understood 



systematic effect that may be related to foreground (see e.g. Ref. (22] ), have been proposed 
(such as e.g. the imprint of the topology of space [H, 24j the breakdown of local isotropy 
due to multiple scalar fields [25[ or the existence of a primordial preferred direction 0, 27|). 

The broken statistical isotropy of the temperature fluctuations may also be related to a 
violation of local isotropy, and thus from a departure from the Friedmann-Lemaitre symme- 
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tries. This can arise either from a late time evolution of the universe (see e.g. Refs. HSSj 
in which it is argued that the subtraction of a Bianchi Vllh leaves a statistically isotropic 
CMB sky) or from the primordial dynamics which would have imprinted the broken statis- 
tical isotropy in the initial conditions. The latter has recently been advocated on the basis 



of a cylindrically symmetric Bianchi / inflationary model [26j. In these models, the shear 
decays as the inverse of the third power of the scale factor so that it can play a significant 
role only in the early stage of the inflationary period. Isotropy is asymptotically reached 
during inflation and the whole subsequent cosmological evolution can be approximated by 
a Friedmann-Lemaitre universe. It follows (1) that the anisotropy is only imprinted in the 
largest wavelengths and (2) that the constraints on the shear of the observable universe from 
the isotropy of the CMB 



B2J, ]33|, |34j or big-bang nucleosynthesis [35| are satisfied. 
The primordial Bianchi / phase modifies the evolution of the modes (in particular gravity 
wave and scalar perturbations shall be coupled through the shear) and initial conditions 
(and thus the quantization procedure) has to be performed in a consistent way with the sym- 
metries of the background spacetime during inflation (see however Ref. (26[ for a proposal 
in a locally rotational invariant and homogeneous spacetime of the Kantowski-Sachs family). 

In this article, we investigate the general theory of gauge invariant perturbations about 
a Bianchi I background spacetime during inflation. Bianchi universes are spatially homoge- 
neous spacetimes and are thus of first importance in cosmology since they express mathe- 
matically the cosmological principle. The study of perturbations in Bianchi / was r oug hed 
out in Ref. [36] where the Bardeen formalism was used (see also Ref. [37| and Ref. 38] for 
the case of higher-dimensional Kaluza-Klein models). A similar work was undertaken in the 



1 + 3 covariant formalism |39j but the identification of gravitational waves and the quantiza- 
tion procedure was not adressed (see Ref. [40(] for the generalisation of the Mukhanov- Sasaki 
variables in this formalism). 

Thus, ^starting from a general parameterisation of the perturbed spacetime a la 
we will define in Section [Til a scalar- vector-tensor decomposition and construct 
Contrary to the Friedmann-Lemaitre case, these three types of 



Bardeen [41 

gauge invariant variables, 
perturbations will be coupled through the shear. In section IIII[ we derive the perturbation 
equations. We then show in Section [TV] that they can be reduced to a set of coupled reduced 
equations with a mixing between scalar and tensor modes; special care will be taken to 
vector modes. This work will allow to generalize the Mukhanov- Sasaki variables and paves 
the way to the study of the cosmological signatures of a primordial anisotropy [49| . 



I. COSMOLOGICAL DYNAMICS OF BIANCHI I UNIVERSES 



Bianchi spacetimes enjoy a group of isometries simply transitive on spacelike hypersur- 



faces (see Refs. [42l . l43l . |44j for a mathematical expositions on Bianchi spacetimes). Thus, 



they are homogeneous. It follows that the cosmic time t is the only essential dynamical coor- 
dinate and Einstein equations will reduce to ordinary differential equations. The dimension 
of their group of isotropy [421 ]. that is the group of isometries leaving a given point fixed, is 
q = 0. 
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A. General form of the metric 



Bianchi I spacetimes are the simplest anisotropic universe models. They allow for differ- 
ent expansion factors in three orthogonal directions. In comoving coordinates, the metric 
takes the general form 

3 

ds 2 = g^dx^dx" = -dt 2 + J2 x ?(t) { dxi Y ■ (1-1) 

i=l 

It includes Friedmann-Lemaitre spacetimes as a subcase when the three scale factors are 
equal. The average scale factor, defined by 

a(t) = [X 1 (t)X 2 (t)X 3 (t)] 1/s , (1.2) 
characterizes the volume expansion. It follows that we can recast the metric (II. ip as 

ds 2 = -dt 2 + a 2 (t)-f ij (t)dx i dx j . (1.3) 

The "spatial metric" 7^ is the metric on constant time hypersurfaces. It can be decomposed 

as 

7ij = exp [2/3i{t)] 5ij , (1.4) 

with the constraints 

3 

X)A = o. (1.5) 

1=1 

Let us emphasize that /3j are not the components of a vector so that they are not subjected 
to the Einstein summation rule. Note also that all latin indices i,j, . . . are lowered with the 
metric 7^. The decomposition (jl.4p implies that 7^ = 2/3,7^-, where a dot refers to a deriva- 
tive with respect to the cosmic time, and it can be checked that the spatial hypersurfaces 
are flat. This relation, together with the constraint fll.5p . implies that the determinant of 
the spatial metric is constant 

7 = = • 

This simply means that any comoving volume remains constant during the expansion of the 
universe, even if this expansion is anisotropic. We define the shear as 



and introduce the scalar b 2 = bija % K This definition is justified from the relation to the 1 + 3 
covariant formalism (see Appendix IA 3h . Let us emphasize at this point that (7^)' = —2a l i 
differs from 7^ = Y p l^ k ipk = +2a t:j . 

Introducing the conformal time as dt = adrj, the metric (11.31) can be recast as 

ds 2 = a 2 (ri) [-drf + 7^(7/) dxW] . (1.7) 

We define the comoving Hubble parameter by TC = a'/ a, where a prime refers to a derivative 
with respect to the conformal time. The shear tensor, now defined as 

1 , . . 
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is related to &ij by cr^- = a&ij. From the relation (7^)' = 7^ = 2/3'^, the definition 

a 2 = a lj( j lj (1.9) 

is explicitely given by 

i=i 

and is related to its cosmic time analogous by a = da. Again, we stress that (j 1 ^) 1 = —2<j^ 
differs from (7')^ = Y p l^ k l' P k = +2o" iJ . 

B. Background equations 

We concentrate on an inflationary phase during which the matter content of the universe 
is assumed to be described by a minimally coupled scalar field, (p, with stress-energy tensor 

T^u = dtfdvtp - (^d a (pd a <p + V) 9fMV . (1.11) 



Making use of the expressions (1A6HA7I) (see Appendix IA 2j) . we easily obtain the Friedmann 
equations 

1 



H =3 



\p' 2 + V{<p)c? 



+ > 2 , (1.12) 




H'=~W 2 -V{ip)a 2 ]-\o 2 (1.13) 

(0$)' = -27toJ , (1.14) 

where k = 8ttG. The first two are similar to the ones usually used in a Friedmann-Lemaitre 
universe, up to the contribution of the shear (which acts as an extra massless field). The 
latter equation arises from the trace- free part of the "if -Einstein equations and gives an 
extra-equation compared to the Friedmann-Lemaitre case. We can easily integrate it and 
conclude that the shear evolves as 

°) = -2 (1-15) 
where <Sj is a constant tensor, («Sj)' = 0. This implies that 

- 2 = ^- 2 = ^, (1.16) 
a 4 a b 

(with S 2 = SjSl) from which we deduce that 

a' = -2Ha. (1.17) 

Let us note that these equations can be combined to give 

2H 2 + W = kcl 2 V , k( v ') 2 = 2H 2 - 2H' - a 2 . (1.18) 

These equations are completed by a Klein-Gordon equation, which keeps its Friedmann- 
Lemaitre form, 

<p" + 2Hp' + a 2 V v = 0. (1.19) 

The general solution for the evolution of the scale factor from these equations is detailed in 
Appendix IA4I 
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II. GAUGE INVARIANT VARIABLES 



This section is devoted to the definition of the gauge invariant variables that describe the 
perturbed spacetime. We follow a method a la Bardeen. In order to define scalar, vector 
and tensor modes, we will need to use a Fourier transform. We start, in § III A[ by recalling 
its definition and stressing its differences with the standard Friedmann-Lemaitre case. In 
§ III B[ we perform a general gauge transformation to identify the gauge invariant variables. 



A. Mode decomposition 

1. Definition of the Fourier transform 

We decompose any quantity in Fourier modes as follows. First, we pick up a comoving 
coordinates system, {x 1 }, on the constant time hypersurfaces. Then, we decompose any 
scalar function as 

f&,v) = I 7—^ / (h,v) e*"* , (2.1) 



with the inverse Fourier transform 



f(k 3 ,v) = J ^j,f(x\ V )e-^\ (2.2) 

In the Fourier space, the comoving wave co-vectors ki are constant, k' { = 0. We now define 
k l = j^kj that is obviously a time- dependent quantity. Contrary to the standard Friedmann- 
Lemaitre case, we must be careful not to trivially identify ki and k l , since this does not 
commute with the time evolution. Note however that Xik % = x l ki remains constant so that 
there is no extra-time dependency entering our definitions fl2.1H2.2p . In the following of 
this article, we will forget the "hat" and use the notation / (x^rj) and / (kj,rj) both for a 
function and its Fourier transform. 

It is easily checked, using the definition (JUBJ), that 

(k*)' = -2a ip k p . (2.3) 

This implies that the modulus of the comoving wave vector, k 2 = k % k^ = ^kikj, is now 
time-dependent and that its rate of change is explicitely given by 



T 

where we have introduced the unit vector 



= -a ij kikj , (2.4) 



k = |. (2.5) 

This vector will turn to be particularly useful for our analysis and we note that it evolves as 

{je)' = (o pq k p k q )k l - 2a ip k p . (2.6) 

Indeed, we find that in the standard Friedmann-Lemaitre limit (<Xy = 0), k % and k are 
constant. 
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2. Decomposition of the vector and tensor modes 

We shall now decompose the perturbations in their scalar, vector and tensor modes. 
Any (3-dimensional) vector field, V 1 , can be decomposed as 

Vi = diV + Vi, with d% = 0, (2.7) 

Note that we have chosen orthogonal (but not Cartesian) coordinates on the (Euclidean) spa- 
tial sections (in particular spatial flatness and homogeneity imply that in these coordinates 
the Christoffel symbols vanish and that dkjij = 0). It follows that its Fourier components 
can be split as 

Vi = kiV + Vi, with k i V i = 0, (2.8) 

so that V 1 lives in the subspace perpendicular to k % . This is a 2-dimensional subspace so that 
Vi has been split into 1 scalar (V) and two vector modes (Vi) that correspond to transverse 
modes. Let us now consider the base {e^e 2 } of the subspace perpendicular to k l . By 
construction, it satisfies the orthonormalisation conditions 

efkjY' = 0, e?e$ 7 ij ' = 5 ab . 

Such a basis is defined up to a rotation about the axis k l . Now, the vector modes can be 
decomposed on this basis as 

Vi(ki,ri)= J2 V a(kv)eKk), (2.9) 

a=l,2 

which defines the two degrees of freedom, V a , which depend on k since the decomposition 
differs for each wave number. The two basis vectors allow to define a projection operator 
onto the subspace perpendicular to k % as 

Pij = e\e) + e\e) = 7o - - hk y (2.10) 

It trivially satisfies P]P 3 k = P%., Pjk^ = and P^jij = 2. It is also the projector on vector 
modes so that we can always make the scalar-vector decomposition 

Vi = >'i}>,-/yi}. (2.ii) 

Analogously, any (3-dimensional) symmetric tensor field, Vy, can be decomposed as 

Vij = T 7lj + Aij S + 2d {i V j} + 2V lj , (2. 12) 

where A^- = didj — A^j/S and 

9^ = 0, V* = = diV ij . (2.13) 

The symmetric tensor is transverse and trace-free. Hence it has only two independent 
components and can be decomposed as 

V l3 (k t , V )= J2 V x (k\ v )e*(k t ) (2.14) 

A=+,x 
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where the polarization tensors have been defined as 



V2 



p^p 2 4- p?p^. 

61+ 1 3 J 1 3 6*. 



V2 



(2.15) 



It can be checked that they are traceless (ehy 1 -* = 0), transverse (e^k 1 = 0), and that the 



two polarizations are perpendicular (e^-ejf = 5*). This defines the two tensor degrees of 
freedom. 

In order to deal with the properties of the polarization tensors, it is useful to define two 
new quantities 

Qij = e\e) - e- e] , and 
The tensor Qy trivially satisfies 



PijQ 



, Q tJ Q 13 = 2 . 

They allow us to simplify the product of two and three polarization tensors as 



(2.16) 
(2.17) 



'ik j 



(P^ + QijT] Xfl ) 







Introducing the projector operator on tensor modes by 



A 



and the "trace extracting" operator 



ab 



pa pb _ p p 



ab 



[2.18) 



T; 



the scalar-vector-tensor terms in the decomposition of Eq. (I2.12p are extracted as follows 



\Vabl 



ab 



Hj + 



V ab T° 



T{j -f- 2k 



(i 



P-)k b V, 



ab 



+ A*°Vab. 



(2.19) 



In this expression, has been split into 2 scalars (T and S), two vector modes (V^) and 
two tensor modes (Vij). Thus, we can always split any equation Vi = by projecting along 
k l (scalar) and P- (vector) and any equation = by projecting along Y 3 (scalar), T JJ 
(scalar), P^ (vector) and A l3 b (tensor). 



3. Properties of the projectors, polarization vectors and tensors 

The previous SVT decomposition matches the one used in the perturbation theory about 
a Friedmann-Lemaitre spacetime. There is however an important difference that we will 
have to deal with. As we pointed out, in a Bianchi I spacetime, the spatial metric is time- 
dependent. It implies in particular that, in order to remain an orthonormal basis perpen- 
dicular to k l during the time evolution, the polarization vectors, and thus the polarization 
tensors, must have a non-vanishing time derivative. Indeed, since (fcj)' = 0, the vector (e*)' 
is orthogonal to ki and is thus a linear combination of e 1 and e 2 , that is 

b 
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In each time hypersurface, there is a remaining freedom in the choice of this basis because 
of the rotational invariance around k l . We can continuously fix the choice of the basis by 
imposing 

Tl[ab] = . 

The orthonormalisation condition implies that (e*e*)' = and thus that 

ft«6 = -<7y44 (2.20) 

Consequently, the time derivative of any polarization vector is given by 

(e?y = J2K ab et + 2a tJ el, (2.21) 

b 

from which we deduce 

k\ef)' = 2a pi k p e1. (2.22) 

This allows us to derive the expression of the time derivative of the polarization tensor. 
Starting from their definitions (12.151) . we easily obtain that 

(4)' = -(^4)^- - (° kl Pki)4 + H4 k > ( 2 - 23 ) 

from which we can deduce some useful algebra 

W(4)' = 0, ''•'( i)' 2rr" i. F(4)' = 2^%4. (2.24) 
We also have that 

(e?)' 4" = . (2.25) 

We gather in Appendix [B] some other useful relations concerning the polarization vectors 
and tensors. 

For the sake of completeness, we shall define here two important matrices for the following 
of our computation, 

l ab = t ij e a e b ) 

which is manifestly symmetric in ab and 



M x ab = 4<e£ , (2.26) 



Kb = QiA4 , (2.27) 

which is anti-symmetric in ah. We stress that a and A are not indices but only labels. It can 
easily be checked that 

^ = 7f(o-°i)^ + 7i(^)^' (2 ' 28) 

and that 



XX- = °- (2-29) 
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B. Defining gauge invariant variables 

1. Gauge invariant variables for the geometry 

Let us consider the most general metric of an almost Bianchi / spacetime. It can always 
be decomposed as 

ds 2 = a 2 [- (1 + 2A) dr] 2 + 2B i dx i dr] + ( 7ii + h {j ) dx'dx*] . (2.30) 

Bi and hij can be further conveniently decomposed as 

Bi = diB + B % , (2.31) 

= 2C fa + ZdifyE + 2d {l E j} + 2E {j , (2.32) 

with 

di& = = diE\ El = = diE ij . (2.33) 

Note that this decomposition of involves the shear. This judicious choice is justified, a 
posteriori, by the simplicity of the transformation properties of the perturbation variables, 
as we shall now see. 

Let us consider an active transformation of the coordinate system defined by a vector 
field £. The coordinates of any point change according to 

x»^3P = x fi - C (x u ) (2.34) 

so that the spacetime metric transforms as 

SV -* 9v» + AfiV > ( 2 - 35 ) 

where C^g^ u is the Lie derivative of g^ u along £. At first order in the perturbations, we 
decompose the metric as g> M „ = + 5g^ v and it follows that 

Sg^u — > $g»v + C^g^. (2.36) 

The vector field £ is now decomposed into a scalar and vector part as 

e = T (x\ V ) , C = d i L(x j ,i]) + V (x\ V ) , (2.37) 

with diL 1 = 0. With the use of the expressions (1C1I) . we deduce that the perturbations of 
the metric transform as (in Fourier space) 

A -> A + T' + HT (2.38) 
(k 2 L)' 

B ^ B -T + ^—r^ (2.39) 
k 

C ^C + HT (2.40) 

E^E + L, (2.41) 



for the scalar variables, and as 



Bi -> Bi + 7ij (U)' - 2ik j a lj P\L (2.42) 
Ei -> + ^ , (2.43) 
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for the vector variables. We also obtain that the tensor modes are readily gauge invariant, 

-> E^ . (2.44) 

Had we not included the shear in the decompostion (12.311) . this would not be the case. 

Let us also note that the transformation rule of the vector modes is different from the 
one derived in Ref. [36] where the non-commutativity between the projection and the time 
evolution has been neglected. 

From the gauge transformations (I2.38H2.4TI) . we can construct a set of gauge invariant 
variables for the scalar sector. Only two degrees of freedom remain, the two other being 
absorbed by the scalar part of the gauge transformation. We define the two gravitational 
potentials 



a 



-C-H 



B 



B 



(k 2 E)' 
' k 2 
(k 2 E)' 
k 2 



(2.45) 



(2.46) 



From the gauge transformations (I2.421l2.43l) . we deduce that a gauge invariant vector per- 
turbation is given by 

$i = - jij (E j )' + 2ik j a lj P l i E, (2.47) 

It is obvious from these expressions that when 7^- is time- independent, that is when 
Oij = 0, these variables reduce to the standard Bardeen variables defined in the Friedmann- 
Lemaitre case. By analogy, we define the Newtonian gauge by the conditions 

B = B l = E = 0, (2.48) 

so that 

A = $, C = -ty, = (2.49) 

the latter condition being equivalent to = —E[ + 2aijEK 



2. Gauge invariant variables for the matter 



We focus our analysis on the scalar field case, which is the most relevant for the study of 
inflation. Under a gauge transformation of the form (12.341) . it transforms as (p — > p> + £%<p. 
At first order in the perturbations, we get 

5<p -> Sip + (2.50) 

that is 

5p — * 5p + = 5<p + <p'T , (2-51) 
with use of Eq. f)2.37p . Thus, we can define the two gauge invariant variables 



Q = S<p- —<?' 



and 



They are related by 



X = 5ip + 



B 



(k 2 E)' 



k 2 



(2.52) 
(2.53) 
(2.54) 
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III. PERTURBATIONS EQUATIONS 



Once the gauge invariant variables have been defined, we can derive their equations of 
evolution. The mode decomposition will require a decomposition of the shear tensor in a 
basis adapted to each wave-number. We start by defining this decomposition and then we 
derive the perturbed Klein-Gordon and Einstein equations. 



A. Decomposition of the shear tensor 

The shear <Jij is a symmetric tracefree tensor and, as such, has 5 independent components. 
In the coordinates system (ll.3til.4p . it was expressed in terms of only two independent 
functions of time A (77). The 3 remaining degrees of freedom are related to the 3 Euler 
angles needed to shift to a general coordinate system. 



1. Components of the shear 

As mentioned before, when working out the perturbations in Fourier space, it would be 
fruitful to decompose the shear in a local basis adapted to the mode we are considering. 
The shear, being a symmetric trace-free tensor, can be decomposed on the basis {ki, e}, e 2 } 
as 

0-y = - ( kkj - -%j ) <7„ + 2 J2 °V<* %e" } + Yl a Tx£ij- (3.1) 

" V 6 J o=l,2 A=+,x 



This decomposition involves 5 independent components of the shear in a basis adapted to 
the wavenumber k{. We must stress however that (a,,,a va ,cr T \) must not be interpreted as 
the Fourier components of the shear, even if they exphcitely depend on fcj. This dependence 
arises from the local anisotropy of space. 

Using Eq. (13. ip . it can be easily worked out that 

<y i3 i ij = o , 

o~ijk % = o~Jij + ^ ~ v a& a j i Oijk % k? — <j|| , 

a 

and 

^ij^X = i 0~ijk €? a = 0~ va . 

The scalar shear is exphcitely given by 

^ = <^ = !< + 2£<. + £^, (3-2) 

"a X 

which is independent of fcj. We emphasize that the local positivity of the energy density of 
matter implies that a 2 /6 < 7i 2 and thus 

i.,<^<«. (3.3) 

This, in turn, implies that 

<r < 2H, (3.4) 
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a property that shall turn to be very useful in the following of our discussion. Analogously, 
we have that 

V6H. (3.5) 

The following derivations will involve the contraction of the shear with the polarization 

vectors, 

(3.6) 
(3.7) 

(3.8) 



A 



A 

ab i 



from which we deduce that 



-a„ . 



It will also involve the contraction of the shear with the polarization tensors, 



A ' 



which implies that 



r, Ji pi _ „ - JiJM _ _ xA 

a il e \^j— cr T\i &U&X j 2 II M ' 



To finish, we will make use of the following expression 



(3.9) 



(3.10) 



2. Time evolution of the components of the shear 

In the previous paragraph we detailed the definition of the components of the shear in a 
basis adapted to the wave mode hi. The time evolution of these modes is easily obtained 
from Eq. (041) 



a 

a' va + 2Ha va = -cr va ^ l{ ~ /J ^vb^xM 



x 

ab i 



b.\ 



a.b 



(3.11) 
(3.12) 
(3.13) 



where the matrix M. x b is defined in Eq. ( 12.261) . 

These equations allow us to derive some important constraints on the rate of change of 
a,, and a T \. Since Eq. (13.111) implies that 



1/2 ^ 

a 2 



2j2< a <° 2 <™ 2 



(3.14) 



we can conclude that 



1 / 2 ^ 



< m. 



(3.15) 
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Identically, Eq. (I3.13P implies that 



a 



-(a 2 a T x)' =2^2M* b a ya <T vb < 



a.b 




so that 



1 



- 2 {a 2 a TX )' <3V2H 2 . 



(3.16) 



a 



The two relations fl3.15H3.16l) will be important at the end of our analysis. 



B. Klein-Gordon equation 




(3.17) 



When expanded at first order in the perturbations, the r.h.s. is trivially given by V v (tp) + 
V<p<p{<p)x- It follows that the Klein-Gordon equation at first order in the perturbations is 
then obtained to be 



where V w is the second derivative of the potential with respect to the scalar field. Surpris- 
ingly, it has the same form as in the Friedmann-Lemaitre case. This can be understood if 
we remind that the d'Alembertian can be expressed as □</? = d v [a/— gg^d^ip] / y/—g, and if 
we realize that \f—g does not involve the shear. Thus, at first order in the perturbations, 
the only place where the shear could appear would be associated with bg % K But then it 
would multiply 8$ which vanishes. Consequently the Klein-Gordon equation is not modi- 
fied. This result is not specific to the scalar field case as the conservation equation in the 
fluid case is also the same as for a Friedmann-Lemaitre spacetime, indeed only as long as 
the anisotropic stress vanishes [see Eq. (1A22j) . 

C. Einstein equations 

The procedure to obtain the mode decomposition of the Einstein equations is somehow 
simple. We start from the general perturbed equation 5G% = k5T^ with the expressions (1C2I - 
IC4I) and (1C10HC121) respectively for the stress-energy tensor and the Einstein tensor and we 
then project them, as described in § III A 21 

Special care must however be taken. In the Friedmann-Lemaitre case, the projections on 
the scalar, vector and tensor modes commute with the time evolution. This no more the 
case in a Bianchi / universe, as explained in § III A3! Let us take an example and consider 
the extraction of the vector part of an equation involving a term of the form ($*)' + 7i$\ 
We project this equation on the polarization tensor e" to get 



We then use Eq. (I2.2ip to rewrite $*(e")', and we develop the shear in the basis adapted to 
the mode k l . This implies, in particular, that contrary to the Friedmann-Lemaitre case, the 
scalar, vector and tensor modes will be coupled. 



X " + 2H X ' - J ij did jX + a 2 V wX = W + 2H<p')$ + <p'(& + 3*') , 



(3.18) 
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This being said, the extraction of the mode decomposition of the Einstein equation is a 
lengthy but straightforward computation that we carry in the Newtonian gauge. It reduces 
to (1) Fourier transforming the Einstein equations, (2) projecting them on the modes, (3) 
commuting the projection operators and the time evolution in order to extract the evolution 
of the polarizations and (4) finally expressing the decomposition of the shear. 



1. Scalar modes 



There are 4 scalar Einstein equations. The first is obtained from SGq = k8Tq and gives 

k 2 Q + 3H{Q' + HQ) - | {^' 2 Q - <p'x' ~ V>\) = 
lk 2 



X = $ + 1) , (3.20) 



\o 2 [X - 3¥] + \-cr^ - \k 2 - \ $>ta£a + «A + 27*^)^,(3.19) 

a A 

where we have defined the extremely useful variable [2} 

Q' 

and the quantity 

a va = ika va . (3.21) 

As an example, the only tricky term which appears when deriving this equation is aj(Ej)', 
which is obtained from 

°){El)' = (a}El)> - Ef(a])' = $> t a£a)' - E\{-2Ho)) = $> tA £ a )' + 2Hct tX E x , 

where we have used Eq. (I1.14p to compute (crj)'. We will not detail these steps in the 
following. 

The second equation is obtained from k l 5G® = nk l 5Tf . We find 

Q' + HQ - ^p'x = "^ 2 * + + \ E a ^ ■ ( 3 - 22 ) 

A 

The two remaining equations are obtained from 

5)5G{ = KS}STi, fk% - 5G{ = k (k% - M 5T? 
and take the form 

+ 2HQ' + HQ' + (2H' + H 2 )Q - \k 2 {Q - tt) + - [p' 2 Q - ip' x ' + V v a 2 x] = 

a 

+ \ Y,K& + «a + 2Ha TX )E x ] , (3.23) 



2 

A 



^ 2 ($ - *) = <7|| 



+ 4fc 2 £ Ai^ b a^a^E x - 2k 2 £ a va $ a . (3.24) 

A, a, ft a 



It can be checked that indeed Eqs. (I3.19[ I3.22ti3.24"]) reduces to their well-known Friedman- 
nian form when the shear vanishes. 
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2. Vector modes 

The two vector equations are obtained from 

e« = 0, heiSG* = . 

They respectively give 

$ a = -2a va X + 4j2 M ab°vbEx, (3.25) 

and 

bX b,X 

+4^JV> v6 (<t t+ 5 a x -a TX 5+)E x , (3.26) 

bX 

where the matrix M a b is defined in Eq. (12.271) . It can be shown that Eq. (13. 26ft results 
from the time derivative of Eq. (13.25!) once Eqs. (13. 11113. 131) are used to express the time 
derivatives of the shear. This a consequence of the Bianchi identities. 



3. Tensor modes 



The equation of evolution of the tensor modes is obtained from e{ 8Gj = 0. To simplify, 
we shall use the shorthand notation (1 — A) for the opposite polarization of A, i.e. it means 
that if A = +, then (1 — A) = x, and vice-versa. With the use of Eq. ( 1B6I) . we obtain 



E'{ + 2HE' X + k 2 E x 



+ 2k 2 Y J M x ab a va ^ h 



a.b 



-2k 2 J2°laEx ~ 2a TX a T+ E (1 _ x) + 2a 2 T{1 _ x) E x . (3.27) 



It can be shown that in the long wavelength limit, the former equations (13.181 13.19H3.241 
I3.25H3.261 13.27ft are equivalent to the ones obtained in a more general gradient expansion of 
Einstein equations on large scales 45 . 



IV. REDUCED EQUATIONS AND MUKHANOV-SASAKI VARIABLES 

The previous equations (I3.18[ I3.19H3.241 13.25H3.261 13. 27|) form a coupled set of equations 
for the scalar, vector and tensor modes. In a Friedmann-Lemaitre spacetime, the three kind 
of perturbations decouple and one can arbitrarily set one of the contributions to zero to 
focus on a given type of mode. This is no more possible here, and in particular, it is not 
possible to neglect the vector modes. Their contribution, as we shall see, is in fact central 
to get the correct set of reduced equations. 

First, we introduce the Mukhanov- Sasaki variables 1(1 for scalar and tensor modes as 



v = aQ , v^/A = aE x , (4.1) 
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exactly in the same way as in a Friedmann-Lemaitre spacetime. These three variables were 
shown to be the canonical degrees of freedom that shall be quantized during inflation when 
a Friedmann-Lemaitre universe is assumed [5fl. 



A. Scalar modes 

Let us introduce these variables in our analysis and start by focusing on the scalar modes. 
First, we note that Eq. (I3.22p can be recast under the more compact form 



(2H - <t„ )X = -if'v + °txE x . (4.2) 



If we now combine Eq. (I3.19P with Eq. (I3.23p . replace the vector mode by its expression (13 . 25[) 
and use the background equations (11.181) . we obtain 

k 2 2 k 2 

HX' + 2(H' + 2H 2 )X + naV^v + k 2 ^ = — ($ - *) + -—a,. . (4.3) 

3 on. 

Now, using Eq. (13.241) to simplify the r.h.s., and again replacing the vector mode by its 
expression (I3.25p . we get 

(m - a }] ) fx' + ^V) + AKa 2 VX + 2mV lf v = Ak 2 ( ^ a 2 Va X - M x ab a Va d Vb E x j . 

^ ' \ a a,b,X J 

(4.4) 

Then, forcing Q in the Klein-Gordon equation (13.181) . using also its background version, we 
obtain 

Q" + 2HQ + k 2 Q + a 2 V w Q + 2a 2 V (p X - <p' (x' + = . (4.5) 

Now, we can replace the last term by using Eq. (14. 4p and the next to last by using Eq. (14 .2p 
to get 



Q" + 2HQ + k 2 Q + a 2 V w Q + 2a 2 V ip X 



(2H - a 



Ak 2 d va x ~ Yl M x ab a Va a Vb E x - Akcl 2 VX - 2/caV't 



a,b,\ 



.(4.6) 



Introducing the definitions (14.11) . we obtain, after some algebra which requires in particular 
Eqs (I3.11H3.131) to express terms such as ^ a b J\A x b &vad'vb, 

This equation is the first central result of this section. 
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B. Tensor modes 



The scalar contribution of the tensor equation (13.271) is exactly given by the relation (14. 4|) . 
so that it reduces, after replacing the vector mode by its expression (I3.25p . to 



2 1 / 2a 2 if'cr T x 



fix = -2/i ( i_A)<x T +cr T x + 2fi X a 2 T{1 _ x) + — I m _ T a I Vkv 



_ \ I iro 



El 2a cr^cr * \ \ i / , . 



a 2 \ 2H — a„ / a 

V \ II / 

This equation is the second central result of this section. 



C. Summary 

We have reduced the perturbation equations to a set of three coupled equations for the 
variables v and fix defined in Eq. (14.1 1) . If we define two new functions z s and z\ by 



Z ' s ^ " a " - aV w + 1 



2a 2 V 2 



a 2 \ 2H- cr„ 



■(^fci) = — + 2al (1 . 
z x a 



1 / 2 X ' 



1 / 2aV, 



2^2 



a 2 I 27Y - <7„ 



(4.9) 



the system reduces to 



v" + { k' 2 - — } v 

z s 



^2 n 2 



1 / 2a?<p'a. 



A* A 



a 2 I 2W-o-„ 

1 /" 2a 2 tp'a T x 
2H-a u 



KV 



+ 



1 / 2a 2 cr TX (T T+ 
2H-a„ 



2cr TX a T _ 



Ai(l-A)- 



Formally, it can be rewritten as 

I/" + fcV + fiV = TV , 
where = (v, //+, /i x ). The matrices f2 and T are defined by 

/ fc 2 - ^ \ 



V" + 



k 2 -^ 
fc 2 - 



/ 




(4.10) 



(4.11) 



(4.12) 



(4.13) 



and the functions K A (?7, ki) and D(r/, fcj) can be read on Eqs. (I4.10H4.1TT) . This is one of the 
central results of our study. 
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When the shear vanishes, these equations decouple and we recover the usual equations [5j 
for the variables v and fi\ so that we only have three physical degrees of freedom. Now, the 
anisotropy of space is at the origin of some interesting effects. First the functions z s and z\ 
are not functions of time only. They depend on ki explicitely through the components of 
the decomposition of the shear. Second, the two types of modes are coupled through a non- 
diagonal mass term. The mass term and the evolution operator cannot be diagonalized at 
the same time so that we expect the equivalent of a see-saw mechanism. The importance of 
the vector modes, that cannot be neglected, has to be emphasized again. Had we neglected 
them, the mass term would not be correct. 



D. Sub-Hubble limit 



Let us consider the behaviour of the mass term appearing in Eqs. (I4.10H4.11~]) in the 
sub-Hubble limit in which k/Ti. 3> 1. We introduce the two slow-roll parameters as 



<P 



6 = 1- 



<p' 2 + 2a 2 V ~ Hep' 

in terms of which the Friedmann equations take the form 

(3 - 6)H<p' + V> 2 = 



— Va 2 + -a 2 
3-e 6 



(4.14) 



(4.15) 



and 



e)H 2 



a' 



We now focus on the behaviour of the functions K^, D, z"/z s and z"/z\ in the sub-Hubble 
regime. We define x = <j/\/&H and use the fact that, since a,. /2 < [see Eq. fl3"4l) ]. 

there exists a < 1 such that < x < a due to the positive energy condition [see Eq. ( 13.31) ] . 
Starting from the definition ( 14. 1 Of) for we have 



— (a 2 a T \) 



2H - cr„ 



+ 2|ct t a| x 



2H - cr„ 



(4.16) 



Now, the property ( 13. 16f) implies that the first term of the right hand side of the inequality 
is smaller than 

2^/k^' 



3V2H 2 



Then, 



2^cp' 



2H - <t„ 



'2e 



2H - tr„ 



H 2 -°i 



n- 



< V2e 



1 + x 
1 — x 



(4.17) 



Now, since x varies in the range < x < a, we deduce that ^/(l + x)/(l — x) < 
a/(1 + «)/(l — a). Eq. (13.51) then implies that the second term of the inequality (14. 161) 
is smaller than 



2VQH 



x 



2H - a„ 
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Then, the absolute value is bounded by 



m - a a 



< 



m 



2H' 
2H 



<7, 



Using the fact that Eq. (J32J) implies |«tJ| < 10H 2 {\a' \ < 6H 2 + \2H<t^\), we obtain that 

2(1 -e) + (1 - f) \/6 + 10" 



< 



y2 yr 



a 



1-5) + 



2W- a 

Gathering all these terms, we thus conclude that 



2(1 - a) 



|Na| < \TeU 



> i / 1 + a 
1 — a 



r(, ^ 2(l-e) + (l-|)y6 + 10' 
6 + 2y3 (l-<5) + - — 1 J 



2(1 -a) 



To summarize, we have shown that 



|N A | < ZW J 



(4.18) 



(4.19) 



where Z is a finite constant. This constant can in principle be quite large since a can be 
arbitrarily close to unity in the worst case of an empty universe. A large Z also corresponds 
to a very ellipsoidal Hubble radius, and this explains why the short wavelength limit has to 
be taken much smaller than the average Hubble radius. 

The same reasoning can be applied for p|, \z"/z B \ and \z"/z\\. Thus, it follows that on 
sub-Hubble scales the three physical degrees of freedom decouple and behave as harmonic 
oscillators, 

V" + k?V = 0. (4.20) 



V. PERTURBATION OF THE ACTION 

In order to construct canonical quantization variables and to properly normalize the am- 
plitude of their quantum fluctuations, one needs to derive the action for the cosmological 
perturbations. We will now demonstrate that the previous equations fl4.10H4.lI1) can be ob- 
tained from the expansion of the action, written in the ADM formalism [46] , at second order. 
Another simpler route would have been to infer the action from the equations of motion, 
which is always possible up to an overall factor, that could then be fixed by considering a 
simple limiting case. Still, we prefer to work out the action at second order since it provides 
a check of the previous computations. 

A. ADM formalism 

In the ADM formalism, we expand the metric as 

ds 2 = — (A 2 — NiN*) dt 2 + 2N i dx i dt + g ij dx i dx j , (5.1) 
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and the Einstein-Hilbert action for a minimally coupled scalar field, takes the form 

S = — I dtd 3 x^/=g [NR® + N (K l:j K ij - K 2 ) - kN (<f J <) ir -<) jr - + 2V (<p)) 
2k J 

+kN~ 1 (<p - JSTdifp)' 4 



(5.2) 



where R^ is the Ricci scalar constructed with the metric g^ and Ka is the extrinsic curva- 



ture, defined as 



Ki 



N~ 



fc-2V ( ,iV i) ), K = K\ 



(5.3) 



Every spatial index is now manipulated with the metric gij. The ADM metric is designed 
in such a way that the constraints arising from the Einstein equations can be immediately 
derived from the action. Varying Eq. (15.21) with respect to the lapse N and the shift Ni, we 
get the Hamiltonian and momentum constraints, respectively 

R {3) - (KijK ij - K 2 ) -2V - Kg ij d t ipdjip + N~ 2 k (ip - N^tf) 2 = , (5.4) 

Vj {K{ - Kdj) - N^k (<p - N j djip) di(p = 0. (5.5) 

Comparing the form ( 12.301) of the metric in Newtonian gauge with Eq. ( 15. ip , we conclude 
that the lapse N and the shift Ni are given by 

iV 2 = (1 + 2$) , N = (5.6) 



and that the metric is 



9ij = a 



Hi - 2^ ( 7« + -£ ) + 2%^) + 2E i3 



(5.7) 



It follows that the Hamiltonian and momentum constraints reduce, at first order, to 



2 1 .. 

-.AM/ - ^d i d j ( - 



3^<x 2 - $ (6H 2 - a 2 - 2 ) 



+-a}di& - a) (E}) - kV v 6<p - K<j>5<p = . 



(5.8) 



and 



* 4- ( — 



+2,% + — AQ^a^djEu-^diEji-KipdiSip = 



(5.9) 

respectively. Once Fourier transformed, written in conformal time and projected along its 
scalar and vector components, we recover precisely Eqs. (I3.19P and (13.221) . 

In order to expand the action up to second order in all first order perturbed quantities, 
we expand the spatial metric as 

9ij = « 2 (lij + hij) . 
The inverse metric and its determinant are then given by 



9 ij = ^ 2 



where 



= -2* ( 7ij - + &ij/H) + 28 (i E 3) + 2E {j 



± + lh + lh 2 -\h)hi 



(5.10) 
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B. Action at zeroth and first orders 



The expansions of the action at zeroth and first orders are 



S 



— [ dtd 3 x [a 3 (-6H 2 + a 2 - 2kV + k^ 2 )] , 
2k J J 



— / dtd 3 x 
2k 



dt v ; 



Si 



— I dtd 3 xa 3 
2k 



i?! 3) + a^kj - 2a ij a i l h jl + + 3* (6H 2 - a 2 



(5.11) 
2kV - Kip 2 ) 



+ $ (6H 2 -a 2 - 2kV - Kip 2 ) - 2kV v 5<p + 2K(p5tp] 



— I dtd 3 x t 0, 



d l (4a*) - <9M ^ 
H 



d 

+ Tt 



A ( -§- ) + a> 3 & ij hij + 12a 3 ff* + 2a 3 k^ 



(5.12) 



where we use the notation X n for the n th order term of the quantity X when expanded in 
perturbations. Note that we have used the background field equations to go from the first 
line to the second line in Eqs. fl5.11ti5.12|) . As can be seen, these two terms can be rewritten 
in terms of total derivatives. It implies that the only nontrivial term will arise from the 
expansion of the action at second order. 

C. Action at second order 



A lengthy but straightforward computation shows that the expansion of the action at 
second order is 



S-2 



1 

2k 



dtd x a 



4 3) + N,R {3) + l-hR {3) + K 2 + \hK x + \h 2 K c 
2 2 8 



1,„ 



--h % MK - NxKt - -N x hK G + Nf}C + « - a^d^Sy - V^Stp 
4 2 V 

-2NtVJ^ - hVJip - hN x V - -h 2 V + -VMV + b(p 2 - 2N x (pb(p 

4 2 J 

+N 2 <p 2 + h<p5<p - hiN^ 2 + h 2 <p 2 - -h)h\^ 2 



(5.13) 
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where 



a 2 R? ] 



4 A 



a ij didj 



2H 



a 2 R. 



(3) 



-d,h lj dM 



-<), I > iJ <)'!>,, - \diha&h li 



VHji 



/C 
/Ci 

/C 2 



-69, + -d% (h jl h fl ) 

-QH 2 + a 2 , 

-2Hh + a ij hij - 2a ij a\h 
2Hhijh ij - IH&ijhPhj - 2a l i h im h m i 



+ 2a ij ajh im l 



+a tJ a lm h im W l - -Ji 2 



(5.14) 

(5.15) 

(5.16) 
(5.17) 

(5.18) 



The construction of the action at second order shall be pursued in Fourier space, since 
many non-local operators appear, such as inverse Laplacian A -1 or (o^SjOjj) -1 , when using 
the constraints. Also, it simplifies the use of the background equations (I3.11H3.13I) for the 
components of the shear a , a va and a T \ which were defined in Fourier space. We recall 
that these components are not the Fourier transforms of the shear but its decomposition in 
a basis adapted to a given mode fcj. 

The integral of any 3-divergence is clearly zero in Fourier space. For instance, let us 
consider a typical term like di(^d 1 ^), then 

J dr 1 d 3 xd l (^d 1 ^) = y"d^d 3 kd 3 q[-k- (k + q)^ k ^ q ]5 (3) (k + q) =0. (5.19) 

Thus, we first express S2 in terms of the Fourier modes and then use conformal time. 
Next, hij is replaced by its expression (15.101) in function of the variables Ej and E^. 
All terms involving Ej either vanish or have the form (E^)', and thus reduce to — $ J in 
Newtonian gauge. Then, we decompose $ J and E^ according to Eqs. (12.91) and (I2.14p . The 
constraint (15. 9p . once expressed in conformal time and in Fourier space, can be projected 
onto its scalar and vector parts in order to obtain the scalar constraint (14.21) and the vector 
constraint (13.251) . Then, we replace <3> a in function of X and E\ using the vector constraint 
and substitute $ by X — \& — (ty/Tl) . We then eliminate X using the scalar constraint. 
Finally, we replace E\ and 5(p by their expressions in terms of the variables fi\ and v [see 
Eq. (JO)]. 

After a tedious calculation, that strictly follows the recipe described just above, the action 
5*2 can be recast under a form that contains only the physical degrees of freedom, 



drjd 3 k{ v'v'* + 



vv 



1 1* 
A* aA* a 



1 



-2<t tX (t t _ 



2a 2 y/K<p'o- TU 



2H 



(v*fJL„ + VHI) 



2a o~ tX o~ t 



2H 




M(1-A)A*A 



(5.20) 



r 



that is, in a more compact way, as 

1 
2 



S-2 



j d7]d 3 k (\v'\ 2 - k 2 \v\ 2 + *v(n - r)v* + r) , (5.21) 
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where T is a total derivative which, for the sake of completeness, is explicitely given in 
Appendix [Dl 

It is clear under this form that the variation of this action with respect to the physical 
degrees of freedom leads directly to the equations of motion (14. 7p and (14.81) . More important, 
it shows that the overall factor is unity. It also follows from this action that the canonical 
momentum associated with v and fi\ are ir v = v'* and tt\ = [i'* x . 



VI. DISCUSSION AND CONCLUSIONS 

In this article, we have presented a full and complete analysis of the theory of cosmological 
perturbations around a homogeneous but anisotropic background spacetime of the Bianchi / 
type. We have described the scalar-vector-tensor decomposition and the construction of 
gauge invariant variables. We have reduced our analysis to a scalar field but it can be easily 
extended to include hydrodynamical matter. 

After presenting the full set of evolution equations for the gauge invariant variables, we 
have shown that the vector modes can be algebraically expressed in terms of scalar and 
tensor modes, so that only three physical degrees of freedom remain, one for the scalar 
sector and two for the tensor sector. Contrary to the Friedmann-Lemaitre case, the scalar, 
vector and tensor perturbation equations do not decorrelate and it was important for the 
consistency of the computation not to neglect the vector modes. We have shown that these 
physical degrees of freedom are the trivial generalization of the Mukhanov-Sasaki variables 
that were derived in a flat Friedmann-Lemaitre universe. 

We have also constructed the action for the cosmological perturbations up to second 
order and demonstrated that, after use of the constraints was made, it only contains the 
physical degrees of freedom and takes a canonical form. We have also shown that in the 
sub-Hubble limit the scalar and tensor degrees of freedom decouple and behave as standard 
harmonic oscillators. It follows that one can apply the standard quantization procedure [B| 
and properly define the normalization of the amplitude of their quantum fluctuations. 

The anisotropy of the underlying space induces two physical effects: (1) the equations 
of motion explicitely involve the wave-number k t and (2) a non-diagonal mass term that 
describes the coupling between scalar perturbation and gravitational waves is at the origin 
of a scalar-tensor see-saw mechanism. 

Since the shear decays as the inverse of the second power of the scale factor, the universe 
isotropizes and tends toward a Friedmann-Lemaitre spacetime. The modes that exit the 
Hubble radius during inflation while the shear is non-negligible will experience the see-saw 
mechanism and will have the primordial anisotropy imprinted on their statistical properties. 
Modes of smaller wavelength will not reflect the anisotropy. It follows that an early Bianchi / 
phase may be at the origin of a primordial anisotrop y o f the cosmological perturbations, 



mainly on large angular scales. The companion article [49| describes such a scenario of early 
anisotropic slow-roll inflation. Since the post-inflationary evolution is well described by a 
Friedmann-Lemaitre spacetime, observable effects, and in particular those related to the 
CMB anomalies we alluded to in the introduction, can be taken into account easily once 
the initial conditions are known. This investigation, that we plan to do later, is beyond the 
scope of the present work. 

Our analysis extends and sheds some light on the robustness of the quantization procedure 
that was developed under the assumption of a Friedmann-Lemaitre background, and thus 
on the predictions of the standard inflationary scenario. We emphasize that this work is 
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very conservative and that no new speculative hypothesis was invoked. Indeed, we are not 
claiming that such a primordial anisotropy is needed. On the one hand, it can be used 
to set stronger constraints on the primordial shear. On the other hand it can also be a 
useful example for the study of second order perturbations, in which a shear appears only 



at first order and induces a correlation between scalar and tensor at second order [47|, |48[ , 
and more generally for the understanding of quantum field theory in curved (cosmological) 
spacetimes |19j . One may for instance wonder whether this analysis can be further extended 
to other Bianchi type or to non-spatially flat spacetimes. 
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APPENDIX A: DETAILS ON BIANCHI I UNIVERSES 



1. Geometrical quantities in conformal time 

Starting from the metric (II. 7p in conformal time, the expressions of the Christoffel sym- 
bols are 

t° 00 = h, r°. H- a ■ a,,. r Qj = m) + a). (Ai) 

where we have used the definition of the shear to express 7^ = la^ so that 

(Y j Y = -%r ij , (A2) 

and indeed trivially (7])' = (5))' = 0. 

We deduce that the non-vanishing components of the Ricci tensor are given by 



a 2 R° = 3H' + a 2 (A3) 
R) = (H' + 2H 2 ) 6) + 2Ha) + {a))', (A4) 



a 2 °* 



where we recall that a 2 = OijO % K The Ricci scalar is 

a 2 R = 6(H' + TL 2 ) + a 2 . (A5) 
The non-vanishing components of the Einstein tensor are thus given by 

a 2 G° = -m 2 + l -a 2 (A6) 

a 2 G) = - \2H' + H 2 + ^a 2 ^j 5) + 2Ha) + (a))'. (A7) 

For a general fluid with stress-energy tensor of the form 

Tpv = pUfiU u + P{g iW + Uf,u u ) + ■Kf.u, (A8) 

where p is the energy density, P the isotropic pressure and tt^ v the anisotropic stress (rr^ u u^ = 
and 7r^ = 0), it implies that the Einstein equation takes the form 

3H 2 = Ka 2 p + ^a 2 , (A9) 

W = -^( P + 3P)-ia 2 , (A10) 
(a))' = -ma) + Ka 2 ^) , (All) 

which correspond respectively to the "00" -component and trace and trace-free part of the 
"^'"-equation. The conservation equation for matter reads 

p' + 3H(p + P) + = , (A12) 

where the r/'-component of tt^ u has been defined as a 2 Tx,ij (so that it) = 7* ft 7ffcj). 

To close this sytem, one needs to specify, as usual, an equation of state for the fluid, that 
is an equation P(p), but also to provide a description for tt^ v . The latter vanishes for a 
perfect fluid and for a scalar field. 
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2. Bianchi I universes in cosmic time 



Starting from the metric (11.31) in cosmic time, the expressions of the Christoffel symbols 



are 



r* - a 1 

Oj 



H8) + ^y*7jy 



The Einstein equations take the form 

3H 2 = 



a 
a 



K P + -a , 
-3Ha] + k%) , 



and the conservation equation for the matter reads 

p + 3H(p + P) + a, - 1 -' = . 



(A13) 

(A14) 

(A15) 
(A16) 

(A17) 



3. Bianchi / universes in the 1 + 3 formalism 

The dynamics of Bianchi universes can be discussed in terms of the 1 + 3 covariant 
formalism (see e.g. Refs. 39, 42| ). This description assumes the existence of a preferred 
congruence of worldlines representing the average motion of matter. The central object is 
the 4-velocity of these worldlines. The symmetries imply that it is orthogonal to the 
hypersurfaces of homogeneity, 

(A18) 



^0- 



The projection operator on the constant time hypersurfaces is defined as 

Its only non-vanishing components being i_y= a 2 {t)^ij{t). 

The central kinematical quantities arise from the decomposition 



+ LO 



(A19) 



For a Bianchi / universe, homogeneity implies that D^f = for all scalar functions (where 
the spatial derivative operator is defined as D^T a =-L„ -L"/ V M 'T a etc.) Since D^P = 0, 
the flow is geodesic and irrotational {ui^ u = 0) so that the acceleration also vanishes, a M = 0, 
and we are just left with the expansion, 0, and the shear E M1 ,. 
It is clear from the form (II. 3p that 



9 = 3H. 



(A20) 



The only non- vanishing components of the shear is expressed simply in terms of the trace-free 
part of the Christoffel symbol as 



V — n z rr 
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so that 

£ 2 = ^ = EA 2 - (A21) 

i=i 

With the general from (1A8[) for the stress-energy tensor, we get the conservation equation 

p + (p + P)0 + S^7T^ = O, (A22) 

which reduces to Eq. (IA12I) . 

The Raychaudhuri [42j equation simplifies to 

6 + -0 2 = -S 2 -4ttG(p + 3P). (A23) 
3 



Since = 3H, the r.h.s. is simply 3d/ a = 3H'/a 2 so that it reduces to Eq. 
The Gauss equation takes the form 



(3) P^ = -u a V a ^ lu - 0S^ + ktt^ + - V ( «P - -0 2 + -S 2 ) . (A24) 



2/11 
3 V ^- 3 02 + 2 : 

We have to be careful here since M Q V Q S M;y is not equal to dtE^. It is given, for the ij- 
component by u a V a Y, i:j = (a 2 ^)' - 2T^ j a 2 cr ik = a 2 [(o"^)' - 2cr ik af\ . 

In the particular case of a Bianchi / spacetime, ^R^ = so that the trace of the 
generalized Friedmann equation 42J reduces to 

Kp _ I 2 + -y? = ®r. 
H 3 2 

Shifting to conformal time, this gives Eq. (1A9j) when ^R = 0. The trace-free part leads 
to Eq. flAlin . Note that this implies that when the anisotropic stress vanishes, a 3 T,^ u is 
constant for the « Q V Q time derivative but that it implies that a 2 <r] is constant in terms of 
the ordinary conformal time derivative. The identification of w^V^ and d t holds only for 
scalars (see e.g. Ref. [471]) . 



4. General solution of the background equations 

It is useful to determine general solutions of the evolution of the background spacetime 0] • 
We concentrate on the particular case in which 7r M „ = (relevant for scalar fields) and first 
set 

A = BiW{t). (A25) 
Equations fll.lOp and fll . 16j) then imply that 

(E B . 2 )^ (i) = 5 or (E*?)^ 2 ^' 

from which we deduce that 




(A26) 
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The constraints (11.51) and (11.101) imply that the E>i must satisfy 



5> = o, J2 B >= s2 



(A27) 



8=1 



i=l 



which are trivially solved by setting 

/2 27T 
-Ssinaj, a p = a + — p, p 6 {1,2, 3}. 

Thus, the general solution is of the form 

m = 

where a is solution of 



(A28) 



jj<S sin (a + ^-i ] x 11". 



3H 2 = k P + 



1&_ 
2 a 6 ' 



(A29) 



(A30) 



Once an equation of state is specified, the conservation equation gives p[a] and we can solve 
for a(t). 

As an example, consider the case of a pure cosmological constant, V = const, and = 0. 
The Friedmann equation takes the form 



H 2 = V n 



a 



with V = kV/3 and a* = (1S/6V0) 1 / 6 . It can be integrated easily to get 



a(t) = a* 



sinh ( 3\/Vnt 



1/3 



(A31) 



Asymptotically, it behaves as the scale factor of a de Sitter universe, a oc exp(v / Vo£) but at 
early time the shear dominates and a oc t 1 ^ 3 . 



APPENDIX B: PROPERTIES OF THE POLARIZATIONS 



We summarize here the main properties of the polarization tensors, defined in § III Al 
The time derivative of the polarization tensor is given by 



-j)k ' 



(4)' = -(^4)^i-(^)4+ 4 4 £ i 

or equivalently, 

( AX\' I kl _A \ jyi I kl tj \ AX , o k A\ 

\ £ j ) = ~\ a £ ki) p j ~ \ a p ki)£j +2<Tj£ k ■ 
In terms of the decomposition (13.11) of the shear tensor, it takes the forms 



:bi) 



(B2) 



-A y 



-a TA Py + ct || 4 + 4ct?.£ A 



(B3) 
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and 

(e?)' = -a TX P; + ^ef + 2a*e£. (B4) 

and the time derivative of the projection operator P i3 is given by 

(jwy = -2a T+ 4 + 2a ^ . (B5) 

Let us also give the following relations that turn to be useful for the derivation of the 
evolution equation of the tensor modes 



lj pi 
a U £ \ a jP £ (l-\) — cr T x°T + 



a 

a a e l la^e%- X) = 0. (B6) 



APPENDIX C: PERTURBED QUANTITIES 

For the sake of completeness, let us give the expression of the Lie derivative (12.361) of the 
displacement £ (12.371) 

Cfg m = -2a 2 (T' + HT) 
C&oi = « 2 K - diT - 2a ^) 

C&ij = a 2 [28^) + 2HT ll3 + 2T(j i3 ] . (CI) 

The expression of the components of the stress-energy tensor of the scalar field at first 
order is 

a 2 5T ° = ifQ-M-V^x, (C2) 
a 2 5T? = -di Wx] , (C3) 
^5Ti = -8) - <p' X ' + Vy X ] ■ (C4) 

Note that they are exactly the same expressions than in a Friedmann-Lemaitre spacetime. 
It comes from the fact that Sg^ never appears. Indeed, the 5T i3 etc. components will be 
different compared to the Friedmann-Lemaitre case. 

In Newtonian gauge the Christoffel symbols at first order take the form 

Ko = & (G5) 
61% = d& (C6) 
1 

— / 

2 



STl = Hh l3 + -h' l3 - 2H$7ij ~ (C7) 



5F oj = hi] - a kj h ki (C8) 

= \{h))' -a kj h ki + h kj a ki 
ST) k = (d 3 h lk + dkh jt - d t h jk ) . (C9) 
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In Newtonian gauge, the expressions of the components of the Einstein tensor at first order 
are 



a 2 5G° n 



a 2 5G] = 5\ 



-2A^ + 6W + 2a 2 m 



a' 



— ) a 2 + —dM 
HJ H 1 

-a^d^j + (Effaj + {6H 2 - a 2 )$ , 



-a 



n 



+ <4d j 



* + * + ( - 



- 2di(^' + H$) 



+-A$, - 2a jk d,E ik + a jk diE jk 



(CIO) 



(Gil) 



2^" + (2H 2 + m') $ + A ($ - *) + 2H& + 



^d j {*-*)-^4d j) d»*+o* 



(C12) 



APPENDIX D: DETAILS CONCERNING THE EXPANSION OF THE ACTION 



The total time derivative T that appears in Eq. (I5.20I) is explicitely given by 

(ip'v) 2 2aip'va ij E ij 



T 



l a\h ml h im + a%jh v H + H(aSipY - 



2H - or,, 



2H - a n 



,_L_^L _ W „ WV - 4aV^» 7 - + ^ U I ^ 



-18^ 2 fta 2 + — + - ' 



n 



n 2 



n 



(Dl) 
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